The parton bubble model (PBM) is connected to the glasma flux tube model 
(GFTM), and predicts the ridge and strong CP violation. 
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In an earlier paper we developed a Parton Bubble Model (PBM) for RHIC, high-energy heavy-ion 
collisions. PBM was based on a substructure of a ring of localized bubbles (gluonic hot spots) which 
initially contain 3-4 partons composed of almost entirely gluons. The bubble ring was perpendicular 
to the collider beam direction, centered on the beam, at midrapidity, and located on the expanding 
fireball surface of Au -|- Au central collisions (0-10%) at y^ijvjv=200 GeV. The bubbles emitted 
correlated particles at kinetic freezeout, leading to a lumpy fireball surface. For a selection of charged 
particles (0.8 GeV/c < pt < 4.0 GeV/c), the PBM reasonably quantitatively (within a few percent) 
explained high precision RHIC experimental correlation analyses in a manner which was consistent 
with the small observed HBT source size in this transverse momentum range. We demonstrated 
that surface emission from a distributed set of surface sources (as in the PBM) was necessary to 
obtain this consistency. A Glasma Flux Tube Model (GFTM) that formed longitudinal flux tubes in 
the transverse plane of two colliding sheets of Color Glass Condensate (CGC), which pass through 
one another had been developed. These sheets create boost invariant flux tubes of longitudinal color 
electric and magnetic fields. A blast wave gives the tubes near the surface transverse flow in the 
same way it gave transverse flow to the bubbles in the PBM. In this paper we consider the equivalent 
characteristics of the PBM and GFTM and connect the two models. The PBM results discussed 
above were all obtained without a jet trigger. Therefore while investigating the ridge in this paper 
we add a jet trigger to the PBM. When one considers a 3-4 GeV/c transverse momentum tagged 
charged particle in combination with another intermediate transverse momentum charged particle 
in the PBM, the ridge correlation is generated and explained. In the GFTM the longitudinal color 
electric and magnetic fields have a non-zero topological charge density FF. These fields cause a local 
strong CP violation which effects charged particle production coming from quarks and anti-quarks 
created in the tube or bubble. For investigating these phenomena we use the PBM alone, without 
a jet trigger. We developed four charged particle correlations which show this strong CP violation 
effect and accumulate from bubble to bubble independent of whether particles are pushed or pulled 
(by the color electric field), and rotated in a right or left handed direction (by the color magnetic 
field). We also show that from previously published analyses of experimental data there is already 
strong evidence for the longitudinal color electric fields, and predict correlations which can be used 
to search for color magnetic field effects. 

PACS numbers: 25.75.Nq, ll.30.Er, 25.75.Gz, 12.38.Mh 



I. INTRODUCTION AND REVIEW OF 
MODELS 



In this paper we discuss the characteristics of the 
parton bubble model (PBM)yy] and connect them 
to related characteristics of the glasma flux tube 
model (GFTM) [3]. We will then proceed to use 
the PBM to discuss and explain the observed ridge 
phenomenon by adding appropriate trigger require- 
ments to the PBM. A successful comparison with 
experimental data is made. We then show how the 
PBM can be used to search for strong CP violations 
or Chern-Simons topological charge We develop 
four charged particle correlations that can be used to 
detect the presence of these color electric and color 
magnetic fields. Then we show that existing pub- 
lished analyses can be used to already provide strong 
experimental evidence for the longitudinal color elec- 



tric fields. Our predicted specific four charged par- 
ticle correlations can be used to search for experi- 
mental evidence for the color magnetic fields. 



A. Parton bubble model development 

Our interest in and the eventual development of 
the PBM goes back many years to the early nineteen 
eighties. Van Hove's work[4] in the early eighties 
considered the bubbles as small droplets of quark- 
gluon plasma (QGP) assumed to be produced in 
ultra-relativistic collisions of hadrons and/or nu- 
clei. His work recognized that experimental detec- 
tors only see what is emitted in the final state at 
kinetic freezeout. Therefore final state predictions 
of his string model work on plasma bubbles or any 
other theory required specific convincing predictions 
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that are observable in the final state which can be 
measured by experiment. His work predicted that 
the final state rapidity distribution dn/dy of hadrons 
would exhibit isolated maxima of width Ay ~ 1 (sin- 
gle bubble) or rapidity bumps of a few units (due to 
those cases producing a few bubbles). These rapid- 
ity regions would also contain traditional signals of 
QGP formation. Even this early work by Van Hove 
needed his bubbles near the surface of the final state 
fireball in order to experimentally measure the par- 
ticles emitted in the final state at kinetic freezeout. 
We and other searched for these Van Hove bubbles, 
but no one ever found any significant evidence for 
them. 

Our next paper Q was our final attempt to make a 
theoretical treatment of the single bubble case (sim- 
ilar to the Van Hove case) which could be experi- 
mentally verified. It is possible that with enough 
statistics one in principle could find single to a few 
bubble events. We developed a number of event gen- 
erators which could possibly provide evidence for 
striking signals resulting from these bubbles. One 
should note that in that paper we had included in 
Sec. 11 mathematical expressions describing how 
charged pions are effected by strong color electric 
and color magnetic fields which are present and par- 
allel in a CP-odd bubble of metastable vacuum. This 
section was motivated by the work of Kharzeev and 
PisarskiQ. For example eq. 3-4 of the section show 
the boosts of the 7r+ and the tt" we estimated. Thus 
in 2000 we were already interested in and were in- 
vestigating and publishing predictions for this phe- 
nomenon. 

When we considered the RHIC quantum interfer- 
ence data in 2000 Jj it became clear that the fire- 
ball surface was rapidly moving outward. This im- 
plied a very large phase space region covered by the 
fireball, and also implied that it would be very un- 
likely for there to be only one isolated bubble (small 
sourcesize) with a large amount of energy sitting on 
the surface. It seemed more likely that there would 
be many bubbles or gluonic hot spots around the 
expanding surface. This led to a paper 8] which 
was our earlier version of Rcf. 1]. In that paper 
we concluded that the behavior of the Hanbury- 
Brown and Twiss (HBT) measurements [7, 9] should 
be interpreted as evidence of a substructure of bub- 
bles located on the surface of the final state fireball 
in the central rapidity region at kinetic freezeout. 
The HBT radii were decreasing almost linearly with 
transverse momentum and implying to us a source 
size of ^2 fm radii bubbles were on the surface and 
could be selected for if we considered transverse mo- 
menta above 0.8 GeV/c. These momenta would al- 
low sufficient resolution to resolve individual bubbles 
of ~2 fm radii. We further concluded that these 



HBT quantum interference observations were likely 
due to phase space focusing of the bubbles pushed 
by the expanding fireball. Thus the HBT measure- 
ments of source size extrapolated to transverse mo- 
menta above 0.8 GeV/c were images of the bubbles. 
The HBT correlation has the property of focusing 
these images on top of each other for a ring of bub- 
bles transverse to and centered on the beam forming 
an average HBT radius. 

Thus our model was a ring of bubbles sitting on 
the freezeout surface with average size ~2 fm ra- 
dius perpendicular to the beam at mid-rapidity. The 
phase space focusing would also lead to angular cor- 
relations between particles emitted by the bubbles 
and be observable. Fig. 1 of that paper shows such 
a bubble geometry. For the background particles 
that account for particles in addition to the bub- 
ble particles we used unquenched HIJING0 with 
its jets removed. We assumed that most of the jets 
were eliminated from the central events because of 
the strong jet quenching observed at RHIC[ll|. We 
investigated the feasibility of using charged parti- 
cle pair correlations as a function of angles which 
should be observable due to the phase space focus- 
ing of the particles coming from individual bubbles. 
At that point we fully expected that these correla- 
tions would be observable in the STAR detector 
at RHIC if one analyzed central Au + Au collisions. 
Correlation analyses are powerful tools in detect- 
ing substructures. Historically substructures have 
played an important role in advancing our under- 
standing of strong (non perturbative) interactions. 
Our earlier paper explained the general character- 
istics of the angular correlation data, and was also 
consistent with HBT measurements in a qualitative 
manner. This motivated us to develop a reason- 
ably quantitative model, the parton bubble modelfH 
which is discussed in the following. 



B. Parton bubble model flj 

In this publication [l| we developed a QCD in- 
spired parton bubble model (PBM) for central (im- 
pact parameter near zero) high energy heavy ion col- 
lisions at RHIC. The PBM is based on a substruc- 
ture consisting of a single ring of a dozen adjoining 
2-fm-radius bubbles (gluonic hot spots) transverse to 
the collider beam direction, centered on the beam, 
and located at or near mid-rapidity. The ring re- 
sides on the fireball blast wave surface (see Fig. 1 of 
Ref.Q). We assumed these bubbles are likely the fi- 
nal state result of quark-gluon-plasma (QGP) forma- 
tion since the energy densities produced experimen- 
tally are greater than those estimated as necessary 
for formation of a quark-gluon-plasma. Thus this is 
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the geometry for the final state kinetic freezeout of 
the QGP bubbles on the surface of the expanding 
fireball treated in a blast wave model. 

The twelve bubble ring creates the average behav- 
ior of bubble formation driven by the energy density 
near the surface of the expanding fire ball of the 
blast wave; that forms the final state surface bubbles 
that emit the final state particles at kinetic freeze- 
out. One should note that the blast wave surface is 
moving at its maximum velocity at freezeout (3c/4). 
For central events each of the twelve bubbles have 
3-4 partons per bubble each at a fixed (j) for a given 
bubble. The transverse momentum (pt) distribution 
of the charged particles is similar to pQCD but has 
a suppression at high pt like the data. 

The bubble ring radius of our model was estimated 
by blast wave, HBT and other general considerations 
to be approximately 8 fm. The bubbles emit corre- 
lated charged particles at final-state kinetic freeze- 
out where we select a pt range (0.8 GeV/c < pt < 
4.0 GeV/c) in order to increase signal to background. 
The 0.8 GeV/c pt cut increases the resolution to al- 
low resolving individual bubbles which have a radius 
of ~2 fm. This space momentum correlation of the 
blast wave provides us with a strong angular corre- 
lation signal. PYTHIA fragmentation functions[13] 
were used for the bubbles fragmentation that gener- 
ate the final state particles emitted from the bubbles. 
The PBM explained the high precision Au -I- Au cen- 
tral (0-10%) collisions at = 200 GeV[li] (the 
highest RHIC energy). The PBM fit to the angu- 
lar correlation data was reasonably quantitative to 
within a few percent (see Sec. 4 of Ref. [l|). 

The correlation functions we employed (like the 
HBT correlation functions) have the property that 
for the difference in angles (difference of momentum 
for HBT) these correlations will image all 12 bub- 
bles on top of each other. This leads to average 
observed angles of approximately 30° in A0, 70° in 
Aj], originating from a source size of about ~2 fm 
radius which is consistent with the HBT correlation. 
Thus the PBM generates A(j) vs A77 charged-particle- 
pair correlations for charged particles with pt in the 
range 0.8 Gev/c to 4.0 GeV/c as displayed in Fig. 
1. See Sec. 4 of Ref. [l| for a reasonably quantita- 
tive successful comparison with data. Furthermore 
the model results were consistent with the Hanbury- 
Brown and Twiss (HBT) observations that the 
observed source radii determined by quantum sta- 
tistical interference were reducing by a considerable 
factor with increasing transverse momentum (pt)- 
The HBT radii are interpreted to reduce from ~6 
fm at Pt ~0.2 GeV/c to -2 fm at pt -1 GeV/c for 
our Pt range. The generally accepted explanation 
for this behavior is that as pt increases radial flow 
increasingly focused the viewed region of the final 



state into smaller volumes. If just one small HBT 
size bubble were emitting all the correlated parti- 
cles, this phenomena would lead to large spikes of 
particles emitted at one limited angular region in 
individual events. This is definitely not observed in 
the Au -I- Au or other collision data at RHIC. There- 
fore, a distributed ring of small sources around the 
beam as assumed in the PBM is necessary to ex- 
plain both the HBT results [§| and the correlation 
data 1J| . The particles emitted from the same bub- 
ble are virtually uncorrelated to particles emitted 
from any other bubble except for momentum con- 
servation requirements. An away side peak in the 
total correlation is built up from momentum conser- 
vation between the bubbles. 

The PBM was also recently extended to 
PBMEil5i] which is identical to the PBM for cen- 
tral collisions (0-5%). For centralities running from 
30-80% jet quenching is not strong enough to make 
jets negligible therefore, a jet component was added 
which was based on HIJING calculations. This jet 
component accounts for more of the correlation as 
one moves toward peripheral bins, and explains all of 
it for the most peripheral collisions. The PBME ex- 
plained in a reasonably quantitative manner (within 
a few percent) the behavior of the recent quantita- 
tive experimental analysis of charge pair correlations 
as a function of centralitvflGj. This further strength- 
ened the substantial evidence for bubble substruc- 
ture. The agreement of the PBM and the PBME 
surface emission models with experimental analyses 
strongly implied that at kinetic freezeout the fire- 
ball was dense and opaque in the central region and 
most centralities (except the peripheral region) in 
the intermediate transverse momentum region (0.8 
< Pt < 4.0 GeV/c). Thus we conclude that the ob- 
served correlated particles are both formed and emit- 
ted from or near the surface of the fireball. In the 
peripheral region the path to the surface is always 
small. 

In the PBM and the PBME the bubbles are pro- 
duced in large numbers, are most symmetric, and 
have unusual features in the most central collisions 
that can be related to the GFTM. Therefore the cen- 
tral region PBM is the model most suitable for com- 
paring to the GFTM. 

The PBM fit to the central production Au + Au 
correlation data[l3| which was observed in the STAR 
TPC detector reasonably quantitatively fit the data. 
The average observed angles were approximately 30° 
in A(j) and 70° in A77, originating from a source 
size of ^2 fm radius. Thus the PBM generates 
A(j) vs At] charged-particle-pair correlations for the 
charged particles with pt in the range 0.8 GeV/c 
to 4.0 GeV/c as displayed in Fig. 1 of the present 
paper. See Sec. 4 of Ref.[l] for a reasonably quanti- 
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tative successful comparison with data (within a few 
percent of the observed correlations). 

The CI correlation used in Fig. 1 and the PBM[1] 
is defined below. The two different types of charged- 
particle-pair correlations are the unlike-sign charge 
pairs (US) and the like-sign charge pairs (LS). The 
charge-independent correlation (CI) is defined as the 
correlation made up of charged-particle-pairs inde- 
pendent of what the sign is, which is the average of 
the US plus the LS correlations. Thus CI = (US 
-I- LS)/2. Therefore the CI is the total charge pair 
correlation observed in the experimental detection 
system within its acceptance. 

We utilize a two particle correlation function in 
the two dimensional (2-D) space^ of versus Ary 
(see Ref.fH Sec. 2.1). The PBM generated the two 
dimensional differences for a pair of particles A0 ver- 
sus A?7 contributing to the correlation function. The 
sum over all pairs coming from the same event and 
then summed over all events is our numerator in the 
correlation function. We then divide this numer- 
ator by a sum over all pairs coming from created 
mixed pairs where, the two particles are produced 
in different events, and then summed over all mixed 
events. The sum over all the particle pair entries 
in the mixed event denominator is normalized to 
equal the sum over all particle pairs in the numera- 
tor. This rescales the CI such that the ratio of the 
numerator divided by the denominator has a mean 
of 1. 

The CI was generated for ^/sNN — 200 GeV cen- 
tral (0-5%) Au -I- Au collisions in the pt range 
0.8 GeV/c to 4.0 GeV/c at RHIC. When the CI 
was corrmared to the corresponding experimental 
analyses a reasonably quantitative agreement 

within a few percent of the observed CI was attained. 
See Ref.[H Sec. 4.2, Ref.^ Sec. V, and Ref.[l6| 
Sec. VI B^. The CI correlation for the 0-5% central- 
ity bin is shown in Fig. 1. This analysis and our pre- 
vious work[l|, [l^ was done without using a jet trig- 
ger. In fact jet production in the central region was 
negligible due to strong jet q uenching pll HTlllsj . As 
one can see in Fig. 1 the central CI correlation of 
A(j) on the near side (At/) < 90°), is sharp and ap- 
proximately jetlike at all A77. In contrast the A?/ de- 
pendence is relatively flat. In the central production 



^ A<j) = <j)i — 'f'2 where <j) is the azimuthal angle of a particle 
measured in a clockwise direction about the beam. Arj = 
Vl ~ V2 which is the difference of the psuedorapidity of the 
pair of particles 

The STAR collaboration preferred to use the conventional 
definition of CI = US + LS which is larger by a factor of 
2 than the CI definition used in this present paper that is 
directly physically meaningful. 



the average observed angles are approximately 30° 
in A(j> and 70° in A77. It is of interest to note that the 
sharp jetlike collimation in A(j) persists at all central- 
ities (0-80%). The elongation in At/ persists in the 
0-30% centrality range and then decreases with de- 
creasing centrality becoming jetlike in the peripheral 
bins. These characteristics are in reasonable quan- 
titative agreement with the fits of the PBM and the 
PBME to experimental data analyses [1 [13, [H, d] ■ 
The difference of the US and LS correlations is de- 
fined as the charge-dependent (CD) correlation (CD 
— US - LS). The 2-D experimental CD correlation 
for centralities (0-80%) has a jetlike shape which is 
consistent with PYTHIA jets (vacuum fragmenta- 
tion) for all centralities. This clearly implies that 
both particle hadroniztion and emission occur from 
the fireball surface region as explained in Ref.[lB] 
Sec. II. 



C. Glasma flux tube model 

A glasma flux tube model (GFTM)[3| that had 
been developed considers that the wavefunctions of 
the incoming projectiles, form sheets of color glass 
condensates (CGC)|l^ that at high energies coUide, 
interact, and evolve into high intensity color elec- 
tric and magnetic fields. This collection of primor- 
dial fields is the Glasma[2^ [2l[, and initially it is 
composed of only rapidity independent longitudinal 
color electric and magnetic fields. An essential fea- 
ture of the Glasma is that the fields are localized 
in the transverse space of the collision zone with a 
size of l/Qs. Qs is the saturation momentum of par- 
tons in the nuclear wavefunction. These longitudinal 
color electric and magnetic fields generate topologi- 
cal Chern-Simons charge|3] which becomes a source 
for particle production. 

The transverse space is filled with flux tubes of 
large longitudinal extent but small transverse size 
^Qj^- Particle production from a flux tube is a 
Poisson process, since the flux tube is a coherent 
state. As the partons emitted from these flux tubes 
locally equilibrate, transverse flow builds due to the 
radial flow of the blast wave[22| . The flux tubes that 
are near the surface of the flreball get the largest 
radial flow and are emitted from the surface. As in 
the parton bubble model these partons shower and 
the higher pt particles escape the surface and do not 
interact. These flux tubes are in the present paper 
considered strongly connected to the bubbles of the 
PBM. The method used to connect the PBM and 
the GFTM is described and discussed in the next 
subsection ID. 

Qs is around 1 GeV/c thus the transverse size of 
the flux tube is about 1/4 fm. The flux tubes near 
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No Trigger 




FIG. 1: The CI (sum of all charged-particle-pairs) correlation for the 0-5% centrality bin with a charged particle pt 
selection (0.8 < pt < 4.0 GeV/c), generated by the PBM. It is plotted as a two dimensional A(j> vs. A77 perspective 
plot. 



the surface are initially at a radius ~5 fm. The an- 
gle wedge of the flux tube is ~l/20 radians or '--'3°. 
Thus the flux tube initially has a narrow range in 
(f). The large width in the Arj correlation which in 
the PBM depended on the large spread in Ar] of the 
bubble partons results from the independent longitu- 
dinal color electric and magnetic fields that created 
the Glasma flux tubes. How much of these longitudi- 
nal color electric and magnetic fields are still present 
in the surface flux tubes when they have been pushed 
by the blast wave will be a speculation of this paper 
for measuring strong CP violation? 

It has been noted that significant features of 
the PBM that generated final state correlations 
which fit the experimentally observed correlation 
datafll, [3, [1^ are similar to those predicted by 
the GFTMQ]. Therefore in the immediately follow- 
ing subsection we assume a direct connection of the 
PBM and the GTM, give reasons to justify it, and 
then discuss its consequences, predictions and suc- 
cesses. 



D. The connection of the PBM and the GFTM 

The successes of the PBM have strongly implied 
that the final state surface region bubbles of the 
PBM represent a significant substructure. In this 



subsection we show that the characteristics of our 
PBM originally developed to fit the precision STAR 
Au -|- Au correlation data in a manner consistent 
with the HBT data; implies that the bubble sub- 
structure we originally used to fit these previous data 
is closely related to the GFTM. 

A natural way to connect the PBM bubbles to the 
GFTM fiux tubes is to assume that the final state 
at kinetic freezeout of a fiux tube is a PBM final 
state bubble. Thus the initial transverse size of a 
fiux tube ^1/4 fm has expanded to the size of ~2 
fm at kinetic freezeout. With this assumption we 
find consistency with the theoretical expectations of 
the GFTM. We can generate and explain the trig- 
gered ridge phenomenon and data (see Sec. II), thus 
implying the ridge is connected to the bubble sub- 
structure. We can predict and obtain very strong ev- 
idence for the color electric field of the glasma from 
comparing our multi-particle charged particle cor- 
relation predictions, and existing experimental cor- 
relation publications (see Sec. III). We have also 
predicted correlations which can be used to search 
for evidence for the glasma color magnetic field (Sec. 
III). 

Of course an obvious question that arises is that 
since a flux tube is an isolated system does a PBM 
bubble, that we assume is the final state of a 
fiux tube at kinetic freezeout, also act as an iso- 
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lated system when emitting the final state parti- 
cles? The near side correlations signals since the 
original PBM[i| have always come virtually entirely 
from particles emitted from the same bubble. Thus 
each bubble has always acted as an isolated system 
similar to the behavior of a flux tube. 



II. THE RIDGE IS FORMED BY THE 
BUBBLES WHEN A JET TRIGGER IS 
ADDED TO THE PBM 



In heavy ion collisions at RHIC there has been 
observed a phenomenon called the ridge which has 
many different explanations [1 [H H [H, [H, [13, [H, 
[29I [sot ■ The ridge is a long range charged particle 
correlation in A77 (very flat), while the A(/) correla- 
tion is approximately jet-like (a narrow Gaussian). 
There also appears with the ridge a jet-like charged- 
particle-pair correlation which is symmetric in Ary 
and A0 such that the peak on the jet-like correla- 
tion is at A77 = and A0 = 0. The At/) correlation of 
the jet and the ridge are approximately the same and 
smoothly blend into each other. The ridge correla- 
tion is generated when one triggers on an interme- 
diate pt range charged particle and then forms pairs 
between that trigger particle and each of all other in- 
termediate charged particles with a smaller pt down 
to some lower limit. The first case we will study 
in this paper is a trigger charged particle between 
3.0 to 4.0 GeV/c correlated with all other charged 
particles which have a pt between 1.1 GeV/c to 3.0 
GeV/c. 

In this paper we will investigate whether the PBM 
can account for the ridge once we add a jet trigger 
to our PBM generator [l|. However this trigger will 
also select jets which previously could be neglected 
because there was such strong quenching 111, EE [3 
of jets in central collisions. A jet trigger had not 
been used in the PBM comparison to all previous 
data. We use HIJING[i3] merely to determine the 
expected number of jets to add for our added jet trig- 
ger. These jet particles were added to our PBM gen- 
erator. Thus our PBM generator now had HIJING 
generated background particles, bubbles of the PBM 
and added jet particles from HIJING. We have al- 
ready shown that our final state particles come from 
hadrons at or near the fireball surface. We reduce 
the number of jets by 80% which corresponds to the 
estimate that only the parton interactions on or near 
the surface are not quenched away, and thus at ki- 
netic freezeout form and emit hadrons which enter 
the detector. This 80% reduction is consistent with 
single 7r° suppression observed in Ref.|l^. We find 
for the reduced HIJING jets that 4% of the Au + 
Au central events (0-5%) centrality at ^/snn = 200 



have a charged particle with apt between 3.0 and 4.0 
GeV/c with at least one other charged particle with 
its Pt greater than 1.1 GeV/c coming from the same 
jet. The addition of the jets to the PBM generator 
provides the appropriate particles which are picked 
up by the trigger in order to form a narrow Arj cor- 
relation signal at which is also a narrow signal in 
A4> at (Fig. 10). This narrow jet signal is present 
in the data and is what remains of jets after 80% are 
quenched away. 

We then form two-charged-particle correlations 
between one-charged-particle with a pt between 3.0 
to 4.0 GeV/c and another charged-particle whose 
Pt is greater than 1.1 GeV/c. The results of these 
correlations are shown in Fig. 2. Fig. 2 is the CI 
correlation for the 0-5% centrality bin with the just 
above stated pt selections on charge pairs. Since we 
know in our Monte Carlo which particles are emit- 
ted from bubbles and thus form the ridge, we can 
predict the shape of the ridge for the above pt cut 
by plotting only the correlation formed from pairs of 
particles that are emitted by the same bubble (see 
Fig. 3). The charge pair correlations are virtually 
all emitted from the same bubble. Those charge 
pair correlations formed from particles originating 
from different bubbles are small contributors which 
mainly have some effect on the away side correla- 
tion. Thus Fig. 3 is the ridge signal which is the 
piece of the CI correlation for the 0-5% centrality 
of Fig. 2, after removing all other pairs except the 
pairs emitted from the same bubble. 



A. Parton pt correlated vs. random 

A very important aspect of the GFTM is the boost 
that the flux tubes get from the radial flow of the 
blast wave. This boost is the same all along the flux 
tube and only depends on how far away from the 
center axis of the blast wave the flux tube is. In 
the PBM the partons in the bubble received a boost 
in Pt from the radial flow field in the blast wave; 
that after final state fragmentation of the bubbles 
gave the generated particles a consistent pt spectrum 
with the data. Since the boost from the blast wave 
depended on the position of the bubble in the ra- 
dial flow field, there should be a correlation between 
partons pt within a bubble. If one redistributed the 
partons with their pt boosts uniformly among the 
bubbles, the over all results for the generated corre- 
lations without a trigger would be unchanged. One 
should note none of our prior work (i.e. PBM 1] or 
PBMEfTB'l) contained a trigger. We have only added 
a trigger here to treat the "triggered ridge" which 
obviously requires it. In our treatment of the GFTM 
(Sec. HI) we have removed the trigger. 
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Triggered particles 




FIG. 2: The CI correlation for the 0-5% centrality bin that results from requiring one trigger particle pt above 3 
GeV/c and another particle pt above 1.1 GeV/c. It is plotted as a two dimensional A(j!> vs. A77 perspective plot. 



Once we require a trigger demanding higher pt 
particles, we start picking out bubbles which have 
more radial flow (harder particles) . Thus correlated 
particles which pass our pt selection and trigger come 
almost entirely from the harder bubbles with more 
radial flow, while the softer bubbles which were sub- 
ject to less radial flow mainly generate low enough 
Pt particles that become background particles or are 
lost to the correlation analysis due to pt selection. 
Figure 2 shows the result of our trigger and pt se- 
lection. If one redistributed the partons with their 
Pt radial flow boost among all the bubbles, then all 
bubbles become equal and there is no longer soft 
and hard bubbles. With this change we can generate 
bubble events. We find that both the non-triggered 
correlation and the A0 correlation are virtually un- 
changed in the region less than about 120° (see Fig. 
4). It seems that we move from the situation where 
we have a few bubbles with a lot of correlated parti- 
cles above our pt cut to a lot of bubbles that have few 
correlated particles above our cuts. However there is 
a difference in the triggered correlation's away side 
or A(/) near the 180° peak. 

In Ref.[l| we discuss this away side effect and at- 
tributed this peak to momentum conservation. This 
is however not the whole story. The away side peak 
is generated by the geometry of the bubble ring 
which requires that for every bubble there is another 
bubble nearly 180° away that is emitting particles. 



These symmetry requirements on the geometry con- 
tribute more to the away side peak especially for 
the case without a trigger or redistributed random 
partons. 

This geometry effect is what causes the away side 
correlation of elliptic flow in a peripheral heavy ion 
collision. In such a collision there is a higher en- 
ergy density aligned with the reaction plane. On 
one side of the reaction plane there are more par- 
ticles produced because of increased energy density, 
and because of the geometry of the situation. On 
the other side of the reaction plane there are also 
more particles produced for the same reasons. Thus 
there is a correlation between particles that are near 
each other on one side of the reaction plane and a 
correlation between particles that are on both sides 
of the reaction plane. Momentum conservation re- 
quires particle communication by bouncing into each 
other, while geometry has no such communication. 
The right hand does not need to know what the left 
hand is doing. However, if physics symmetry makes 
them do the same thing then they can show a cor- 
relation. 



B. The away side peak 

The away side peak depends on the fragmenta- 
tion of the away side bubble. For the triggered case 
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Ridge signal 




FIG. 3: The ridge signal is the piece of the CI correlation for the 0-5% centrality of Fig. 2 after removing all other 
particle pairs except the pairs that come from the same bubble. It is plotted as a two dimensional vs. Ar; 
perspective plot. As explained in the text charge pair correlations formed by particles emitted from different bubbles 
are small contributors which mainly have some effect on the away side correlation for /S.(j) greater than about 120° . 



where there is a hard bubble with a strong corre- 
lation of parton pt inside the bubble (GFTM like), 
we will have more correlated particles adding to the 
signal on the near side. On the other hand the away 
side bubble will be on average a softer bubble which 
will have a lot less particles passing the pt cuts and 
thus have a smaller signal. If we consider the case 
of a random parton pt bubble the particles passing 
the Pt cuts should be very similar on the near and 
the away side. In Fig. 4 we compare the two cor- 
relations generated by the correlated and the ran- 
dom cases. We compare the A0 correlation for the 
two cases (solid is correlated and open is random)for 
each of the five A?7 bins which cover the entire A77 
range (0.0 to 1.5). The vertical correlation scale is 
not offset and is correct for the largest Ary bin, 1.2 
< Ary < 1.5, which is the lowest bin on the figure. 
As one proceeds upwards to the next A77 bin the 
correlation is offset by -1-0.05. This is added to the 
correlation of each subsequent Ar; bin. The smallest 
A?7 bin on top of Fig. 4 has a -1-0.2 offset. A solid 
straight horizontal line shows the offset for each A?7 
bin. Each solid straight horizontal line is at 1.0 in 
correlation strength. We see that the away side cor- 
relation in the A^ region greater than about 120° is 
larger for the random case. 



C. Predicted transverse momentum 
dependence of correlation 



Particles which come from bubbles that satisfy the 
trigger (3.0 < pt < 4.0 GeV/c) are harder and have a 
different pt distribution compared to the case with- 
out trigger requirements. In a given triggered event 
if one could select particles that come from the bub- 
ble or ridge, one would find that the pt spectrum is 
harder than the average pt spectrum. In Fig. 5 we 
show how one can by trigger choose particles that 
are rich with ridge particles. We consider a typi- 
cal trigger particle and the associated regions in the 
available experimental (/) and rj ranges. The ridge 
particles we consider are separated from the trigger 
particle by 0.7 in rj. This is done to eliminate par- 
ticles coming from jet production that could also be 
associated with the trigger particle. The bulk of the 
ridge particles lie within 20° of the (j) of the trigger 
particle. The width of the spread is At/) = 40°. 
Therefore on a trigger by trigger basis we accumu- 
late particles from the ridge region (see Fig. 5) and 
form a pt spectrum of charged-particles-pairs. We 
then form a ratio of the pt spectrum for the ridge 
region particles divided by the pt spectrum for all 
particles in central Au + Au events. This ratio is 
highly dependent on our ridge cut area. We can vir- 
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FIG. 4: The A<j!> CI correlation for the 0-5% centrality bin requiring one trigger particle pt above 3 GeV/c and another 
particle pt above 1.1 GeV/c. We compare the two correlations generated by the correlated and the random cases. 
We compare the A(j> correlation for the two cases (solid is correlated and open is random) for each of the five Ar; bins 
which cover the entire A-q range (0.0 to 1.5). The vertical correlation scale is not offset and is correct for the largest 
Arj bin, 1.2 < A77 < 1.5, which is the lowest bin on the figure. As one proceeds upwards to the next Arj bin the 
correlation is offset by +0.05. This is added to the correlation of each subsequent Arj bin. The smallest Arj bin on 
top of Fig. 4 has a +0.2 offset. A solid straight horizontal line shows the offset for each Arj bin. Each solid straight 
horizontal line is at 1.0 in correlation strength. We see that the away side correlation for A<j) greater than about 120° 
is larger for the random case. 



tually remove this dependency by normalizing the 
two Pt spectra to have the same counts for the 2.0 
GeV/c bin. Figure 6 is this ratio which can be con- 
sidered a correlation of particles as a function of pt 
in the ridge region compared to particles from the 
event in general. In Fig. 6 we form this ratio or 
correlation for the bubbles which have a correlated 
Pt among the partons (GFTM like). We also show 
the correlation for bubbles which have a random pt 
distribution among partons in the bubble. For the 
random case the pt spectrum does not differ much 
from the average pt spectrum while there is a big dif- 
ference for the flux tube like case. This correlation 



could be easily measured in the RHIC data. 



D. Comparison to data 

Triggered angular correlation data showing the 
ridge was presented at Quark Matter 2006 31| . Fig- 
ure 7 shows the experimental vs. A77 CI corre- 
lation for the 0-10% central Au + Au collisions at 
y/SNN — 200 GeV; requiring one trigger particle pt 
between 3 and 4 GeV/c and an associated particle 
Pt above 2.0 GeV/c. The yield is corrected for the 
finite A77 pair acceptance. For the PBM generator. 
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FIG. 5: In this figure we display the (p and rj ranges considered for producing ridge particles in the An + An central 
collisions. A typical trigger particle {pt between 3.0 to 4.0 GeV/c) implies two associated regions in the considered (j) 
and ?7 ranges where the ridge particles lie. The ridge particles we consider are separated by 0.7 in ri (Arj) so that one 
eliminates particles coming from jet production. The bulk of the ridge particles lie within 20° of the (j> of the trigger 
particle. The width of the <j> spread is A<j> = 40°. 



we then form a two-charged-particle correlation be- 
tween one charged particle with a pt between 3.0 to 
4.0 GeV/c and another charged particle whose pt is 
greater than 2.0 GeV/c. The results of this correla- 
tion is shown in Fig. 8. Figure 7 shows the corrected 
pair yield determined in the central data whereas 
Fig. 8 shows the correlation function generated by 
the PBM which does not depend on the number of 
events analyzed. We can compare the two figures, 
if we realize that the away side ridge has around 
420,000 pairs in Fig. 7 while in Fig. 8 the away 
side ridge has a correlation of around 0.995. If we 
multiply the correlation scale of Fig. 8 by 422,111 
in order to achieve the number of pairs seen in Fig. 
7, the away side ridge would be at 420,000 and the 
peak would be at 465,000. This would make a good 
agreement between the two figures. 



We know in our Monte Carlo which particles come 
from bubbles and which particles form the ridge. 
The correlation formed by the ridge particles is gen- 
erated almost entirely by particles emitted by the 
same bubble. We have shown in all our publica- 
tions that the same side correlation signals are al- 
most entirely formed by particles coming from the 
same bubble. Thus we can predict the shape and 
the yield of the ridge for the above pt trigger selec- 
tion and lower cut, by plotting only the correlation 
coming from pairs of particles that are emitted by 
the same bubble (see Fig. 9). 

In Ref. 31] it was assumed that the ridge yield 
was flat across the acceptance while in Fig. 9 we see 
that this is not the case. Therefore our ridge yield is 
approximately 35% larger than estimated in Ref. [3l| . 
Finally we can plot the jet yield that we had put into 
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FIG. 6: On a trigger by trigger basis we accumulate particles from the ridge region (see Fig. 5). We then form a 
ratio the of pt spectrum for the ridge region particles divided by the pt spectrum for all particles in central Au + 
Au events. We normalize the two pt spectra to have the same counts in the 2.0 GeV/c bin. This ratio or correlation 
for the bubbles which have a correlated pt among the partons (GFTM like) are the solid points. The random pt 
distribution among partons in the bubble are the open points. 



our Monte Carlo. We used HIJING[T31 to determine 
the number of expected jets, and then reduced the 
number of jets by 80%. This assumes that only the 
parton interactions on or near the surface that form 
hadrons at kinetic freezeout are not quenched away 
and thus enter the detector. This 80% reduction 
is consistent with single tt" suppression observed in 
Ref ■ [isj . This jet yield is plotted in Fig. 10 where 
we subtracted contributions from the bubbles and 
the background particles from Fig. 7. 



III. STRONG CP VIOLATIONS OR 
CHERN-SIMONS TOPOLOGICAL CHARGE 

A. The Source FF. 

The strong CP problem remains one of the most 
outstanding puzzles of the Standard Model. Even 
though several possible solutions have been put for- 
ward it is not clear why CP invariance is respected 
by the strong interaction. It was shown however 
through a theorem by Vafa- Witten [s^l that the true 
ground state of QCD cannot break CP. The part of 
the QCD Lagrangian that breaks CP is related to 
the gluon-gluon interaction term FF. 
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^/spfN = 200 GeV requiring one trigger particle pt between 3 to 4 GeV/c and an associated particle pt above 2.0 
GeV/c. 



Triggered particles 




FIG. 8: The PBM generated CI correlation for the 0-5% centrality bin requiring one trigger particle pt above 3 
GeV/c and another particle pt above 2.0 GeV/c plotted as a two dimensional A(/!> vs. A?; perspective plot. The 
trigger requirements on this figure are the same as those on the experimental data in Fig. 7. 
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Ridge 




FIG. 9: The ridge signal is the piece of the CI correlation for the 0-5% centrality of Fig. 8 after removing all other 
particle pairs except the pairs that come from the same bubble. It is plotted as a two dimensional vs. At; 
perspective plot. 



Jet Signal 




FIG. 10: The jet signal is left in the CI correlation after the contributions from the background and all the bubble 
particles are removed from the 0-5% centrality (with trigger requirements) of Fig. 8. It is plotted as a two dimensional 
A<^ vs. Arj perspective plot. 
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The part of FF related to CP violations can 
be separated into a separate term which then can 
be varied by multiplying this term by a parame- 
ter called 9. For the true QCD ground state 9 = 
(Vafa-Witten theorem). In the vicinity of the 
deconfined QCD vacuum metastable domains 33] 9 
non-zero could exist and not contradict the Vafa- 
Witten theorem since these are not the QCD ground 
state. The metastable domains CP phenomenon 
would manifest itself in specific correlations of pion 
momenta d, [S^. 



B. Pionic measures of CP violation. 

The glasma flux tube model (GFTM) [2] considers 
the wavefunctions of the incoming projectiles, form 
sheets of CGCJJj at high energies that collide, in- 
teract, and evolve into high intensity color electric 
and magnetic fields. This collection of primordial 
fields is the Glasmafl^, HH. Initially the Glasma 
is composed of only rapidity independent longitudi- 
nal (along the beam axis) color electric and mag- 
netic fields. These longitudinal color electric and 
magnetic fields generate topological Chern-Simons 
chargeQ through the FF term and becomes a source 
of CP violation. How much of these longitudinal 
color electric and magnetic fields are still present in 
the surface flux tubes when they have been pushed 
by the blast wave is a speculation of this paper for 
measuring strong CP violation? The color electric 
fleld which points along the fiux tube axis causes an 
up quark to be accelerated in one direction along 
the beam axis, while the anti-up quark is acceler- 
ated in the other direction. So when a pair of quarks 
and anti-quarks are formed they separate along the 
beam axis leading to a separated tt"*" tt" pair along 
this axis. The color magnetic field which also points 
along the fiux tube axis (which is parallel to the 
beam axis) causes an up quark to rotate around the 
flux tube axis in one direction, while the anti-up 
quark is rotated in the other direction. So when a 
pair of quarks and anti-quarks are formed they will 
pickup or lose transverse momentum. These changes 
in pt will be transmitted to the tt"*" tt~ pairs. 

It is important to note that the CP violating asym- 
metries in 7r+ and 7r~ momenta arise through the 
Witten- Wess-Zumino^^, [3^ term. The quarks and 
anti quarks which, later form the 7r+ and it~ , di- 
rectly respond to the color electric and color mag- 
netic fields and receive their boosts at the quark and 
anti quark level before they are pions. These boosts 
are transmitted to the pions at hadronization. Thus 
no external magnetic field is required in the method- 
ology followed in this paper. The following references 
demonstrate this:^^ @, [ssl [36j . 



To represent the color electric field effect we as- 
sume as the first step for the results being shown 
in this paper; that we generate bubbles which have 
an added boost of 100 MeV/c to the quarks in the 
longitudinal momentum which represents the color 
electric effect. The 7r+ and the 7r~ which form a pair 
arc boosted in opposite directions along the beam 
axis. In a given bubble all the boosts are the same 
but vary in direction from bubble to bubble. For 
the color magnetic field effect, we give 100 MeV/c 
boosts to the transverse momentum in an opposite 
way to the tt"*" and the 7r~ which form a pair. For 
each pair on one side of a bubble one of the charged 
pions boost is increased and the other is decreased. 
While on the other side of the bubble the pion for 
which the boost was increased is now decreased and 
the pion which was decreased is now increased by the 
100 MeV/c. All pairs for a given bubble are treated 
in the same way however each bubble is random on 
the sign of the pion which is chosen to be boosted 
on a given side. This addition to our model is used 
in the simulations of the following subsections. 



C. Color electric field pionic measure. 

Above we saw that pairs of positive and nega- 
tive pions should show a charge separation along the 
beam axis due to a boost in longitudinal momen- 
tum caused by the color electric field. A measure of 
this separation should be a difference in the pseudo- 
rapidity (A77) of the opposite sign pairs. This At^ 
measure has a well defined sign since we defined this 
difference measuring from the it~ to the tt"*". In or- 
der to form a correlation we must pick two pairs for 
comparison. The pairs have to come from the same 
bubble (a final state of an expanded flux tube) since 
we have shown by investigating the events generated 
by the model that pairs originating from different 
bubbles will not show this correlation. Therefore we 
require the 7r+ and the tt" differ by 20° or less in (j). 
Let us call the first pair A771. The next pair (A772 
having the same (j) requirement) has to also lie inside 
the same bubble to show this correlation. Thus we 
require that there is only 10° between the average (j) 
of each bubble. This implies that at most in 4> no 
two pions can differ by more than 30°. In Fig. 11 we 
show two pairs which would fall into the above cuts. 
Aryi and A7y2 are positive in Fig. 1 1 . However if we 
would interchange the tt"*" and tt^ on cither pair the 
value of its A77 would change sign. Finally the mean 
value shown on Fig. 11 is the mid-point between the 
TT"*" and the 7r~ where one really uses the vector sum 
of the TT"*" and the 7r~ which moves this point toward 
the harder pion. 

Considering the above cuts we defined a correla- 
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tion function where we combine pairs each having a 
A77. Our variable is related to the sum of the abso- 
lute values of the individual A?7's ( | A771 1 + 1 Ary2 1 ) • We 
assign a sign to this sum such that if the sign of the 
individual A?7's are the same it is a plus sign, while if 
they are different it is a minus sign. For the flux tube 
the color electric field extends over a large pseudora- 
pidity range therefore let us consider the separation 
of pairs | Aryj greater than 0.9. For the numerator of 
the correlation function we consider all combinations 
of unique pairs (sign (lAryij + |A?72|)) from a given 
central Au + Au event divided by a mixed event de- 
nominator created from pairs in different events. We 
determine the rescale of the mixed event denomina- 
tor by considering the number of pairs of pairs for 
the case \A<j>\ lying between 50° and 60° for events 
and mixed events so that the overall ratio of this 
sample numerator to denominator is 1. By picking 
50° < |A0| < 60° we make sure we are not choosing 
pairs from the same bubble. For a simpler notation 
let (sign (|Aryi| -I- |A?72|)) = Ar/i + A?72 which varies 
from -4 to -|-4 since we have an over all rj acceptance 
-1 to +1 (for the STAR TPC detector for which we 
calculated). The value being near ± 4 can happen 
when one has a hard pion with pt of 4 GeV/c (upper 
cut) at r/ = 1 with a soft pion pt of 0.8 GeV/c (lower 
cut) at ?7 = -1 combined with another pair; a hard 
pion with pt of 4 GeV/c at 77 = -1 with a soft pion 
with Pt of 0.8 GeV/c at ?/ = 1. 

In Fig. 12 we show the correlation function of op- 
posite sign charge-particle-pairs paired and binned 
by the variable A771 + A?72 with a cut |A?7| greater 
than 0.9 between the vector sums of the two pairs. 
The events are generated by the PBM[i| and are 
charged particles of 0.8 < pt < 4.0 GeV/c, and |?7| 
< 1, from Au + Au collisions at ^/snn — 200 GeV. 
Since we select pairs of pairs which are near each 
other in cj) they all together pick up the bubble cor- 
relation and thus these pairs of pairs over all show 
about a 4% correlation. In the 1.0 < Ar]i+Ar]2 < 2.0 
region the correlation is 0.5% larger than the -2.0 < 
A771 +Ar]2 < -1.0 region. This means there are more 
pairs of pairs aligned in the same direction compared 
to pairs of pairs not aligned. This alignment is what 
is predicted by the color electric field effect presented 
above. In fact if one has plus minus pairs all aligned 
in the same direction and spread across a pseudo- 
rapidity range, locally at any place in the pseudo- 
rapidity range one would observe an increase of un- 
like sign charge pairs compared to like sign charge 
pairs. At small A77 unlike sign charge pairs are much 
larger than like sign charge pairs in both the PBM 
and the data which agree. See Figs. 10, 11, and 14 
of Ref.[l|. Figure 11 of Ref. 1] compares the total 
correlation for unlike-sign charge pairs and like-sign 
charge pairs in the precision STAR central produc- 



tion experiment for Au + Au central collisions (0- 
10% centrality) at y/JNN=200 GeV, in the trans- 
verse momentum range 0.8 < pt < 2.0 GeV/c[l3|. 
The unlike-sign charge pairs are clearly larger in the 
region near A0 = Arj = 0.0. The increased cor- 
relation of the unlike-sign pairs is 0.8% ± 0.002%. 
Figure 10 of Ref.[H shows that the PBM fit to these 
data gives the same results. Figure 14 of Ref.[l| 
shows that the CD — unlike-sign charge pairs minus 
like-sign charge pairs is positive for the experimental 
analysis. Therefore the unlike-sign charge pairs are 
considerably larger than the like-sign charge pairs. 
In fact this effect is so large and the alignment is so 
great that when one adds the unlike and like sign 
charge pairs correlations together there is still a dip 
at small A77 and A(/) see Fig. 1 and Fig. 3 of the 
present paper. The statistical significance of this 
dip in the two high precision experiments done inde- 
pendently from different data sets gathered 2 years 
apart [3, [3 

is huge. It would require a fluctuation 
of 14cr to remove the dip. This dip is also not due to 
any systematic error since both of the just cited pre- 
cision experiments carefully investigated that possi- 
bility; and found no evidence to challenge the reality 
of this dip. This highly significant dip (~14(t) means 
that like-sign pairs are removed as one approaches 
the region A77 = A0 — 0.0. Thus this is very strong 
evidence for the predicted effect of the color electric 
field. 



D. Color magnetic pionic measure. 

After considering the color electric effect we turn 
to the color magnetic effect which causes up quarks 
to rotate around the flux tube axis in one direction, 
while the anti-up quarks rotate in the other direc- 
tion. So when a pair of quarks and anti-quarks are 
formed they will pickup or lose transverse momen- 
tum. These changes in pt will be transmitted to 
the 7r+ and 7r~ pairs. It was previously shown in 
Sec. Ill A that the CP violating asymmetries in the 
7r+ and the 7r~ momenta arise through the Witten- 
Wess-Zumino term. The quarks and anti quarks 
which, later form the 7r+ and the tt^, directly re- 
spond to the color electric and color magnetic fields 
and receive their boosts at the quark and anti quark 
level before they hadronize into pions. These boosts 
are transmitted to the 7r+ and the 7r~ . 

In order to observe these differential pt changes 
one must select pairs on one side of the bubble in (j) 
and compare to other pairs on the other side of the 
same bubble which would lie around 40° to 48° away 
in (f>. We defined a pair as a plus particle and minus 
particle with an opening angle {0) of 16° or less. We 
are also interested in pairs that are directly on the 
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FIG. 11: Pairs of pairs selected for forming a correlation due to the boost in longitudinal momentum caused by the 
color electric field. The largest separation we allow in (f> for plus minus pair 1 is A0i = 20°. This assures the pair is 
in the same bubble. The same is true for pair 2 so that it will also be contained in this same bubble. The mid-point 
for pair 1 and 2 represents the vector sum of pair 1 and 2 which moves toward the harder paxticle when the momenta 
differ. These mid-points can not be separated by more than 10° in A(j) in order to keep all four particles inside the 
same bubble since the correlation function is almost entirely generated within the same bubble. The Ar/ measure is 
the angle between the vector sum 1 compared to the vector sum 2 along the beam axis (for this case |Ar/| = 1.1). 
The positive sign for Arji comes from the fact that one moves in a positive rj direction from negative to positive. The 
same is true for Ar?2. If we would interchange the charge of the particles of the pairs the sign would change. 



other side of the bubble. We require they are near in 
pseudorapidity (A?7 < 0.2). The above requirements 
constrain the four charged particles comprising both 
pairs to be contained in the same bubble and be close 
to being directly on opposite sides of the bubble (a 
final state expanded flux tube). The difference in pt 
changes due to and predicted by the color magnetic 
effect should give the 7r+ on one side of the bubble 
an increased pt and a decreased pt on the other side, 
while for the 7r~ it will be the other way around. 
This will lead to an ant i- alignment between pairs. 
In Fig. 13 we show two pairs which would fall into 
the above cuts. Both pairs are at the limit of the 



opening angle cut 9i and 02 equal 16°. The pt of the 
plus particle for pair number 1 is 1.14 GeV/c, while 
the minus particle is 1.39 GeV/c. Thus APu is equal 
to -0.25 GcV/c. The pi of the plus particle for pair 
number 2 is 1.31 GeV/c. the minus particle is 0.91 
GeV/c and APta is equal to 0.40 GeV/c. Finally 
the mean value shown on Fig. 13 is the mid-point 
betwcxai the 7r+ and the tt^ where one really uses 
the vector sum of the tt"*" and the tt^ which moves 
this point toward the harder pion. 

Considering the above cuts we defined a correla- 
tion function where we combine pairs each having a 
APt. Our variable is related to the sum of the abso- 
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FIG. 12: Correlation function of pairs formed to exhibit the effects of longitudiual nionientuni boosts by the color 
electric field as defined in the text. This correlation shows that there are more aligned pairs (correlation is larger 
by ^0.5% between 1.0 < Ar?! + Arj2 < 2.0 compared to between -2.0 < Ar]i + Ar]2 < -1.0). Arji + A772 is equal 
to |Ar?i| + |Ar72|. As explained in the text this means there are more pairs of pairs aligned in the same direction 
compared to pairs not aligned as predicted by the color electric field. The signs and more detail axe also explained 
in the text. 



lute values of the individual APj's {\APn\ + 1 APtal). 
We assign a sign to this sum such that if the sign of 
the individual APt's arc the same it is a plus sign, 
while if they are different it is a minus sign. For 
the flux tube the color magnetic field extends over 
a large psciidorapidity range where quarks and anti- 
quarks rotate around the flux tube axis, therefore 
we want to sample pairs at the different sides of the 
tube making a separation in (j) {A(j>) between 40° to 
48°. We are interested in sampling the pairs on the 
other side so we require the separation in 77 (A77) 
be 0.2 or less. For the numerator of the correlation 
function we consider all combinations of unique pairs 
(sign (|APii| + |APj2|)) from a given central Au + 
Au event divided by a mixed event denominator cre- 
ated from pairs in different events. We determine the 



rescale of the mixed event denominator by consider- 
ing the number of pairs of pairs for the case \Ar]\ 
lying between 1.2 and 1.5 plus any value of |A(/)| for 
events and mixed events so that the overall ratio of 
this sample numerator to denominator is 1. By pick- 
ing this Arj bin for all \A4>\ we have aroimd the same 
pair count as the signal cut with the A(j) correlation 
of the bubbles being washed out. For a simpler no- 
tation let (sign (|APn| + |APt2|)) = APti + APta 
which we plot in the range from -4 to -f 4 since we 
have an over all pt range 0.8 to 4.0 GeV/c. Thus the 
maximum magnitude of AP^'s is 3.2 GeV/c which 
makes AP^i -|- APt2 have a range of ±6.4. However 
the larger values near these range limits occur very 
rarely. 

In Fig. 14 we show the correlation function of op- 
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FIG. 13: Pairs of pairs selected for forming a correlation exhibiting the effect of changes in pt due to the color 
magnetic field. The largest opening angle 8 for plus mirms pairs is 16° or less. This opening angle assures that each 
pair has a high probability that it arises from a quark anti-quark pair. In this figure we have picked two pairs at 
this limit {Oi = 16° and Q-i — 16°). The mid-point for pair 1 and 2 represents the vector sum of pair 1 and 2 which 
moves toward the harder particle when the momenta differ. These mid-points are chosen to have 40° < |A0| < 48° 
in order for the pairs to be on opposite sides of the bubble. Since we are interested in pairs directly across the bubble 
we make the At; separation be no more than 0.2. The difference in pt for pair 1 is APti ~ -0.25 GeV/c, while the 
difference in pt for pair 2 is APt2 = 0.40 GeV/c. The minus sign for 1 follows from the fact that the plus particle has 
1.14 GeV/c and the minus particle has 1.39 GeV/c. The plus sign for 2 follows from the fact that the plus particle 
has 1.31 GeV/c and the minus particle has 0.91 GeV/c. 



posite sign charged-particle-pairs paired and binned 
by the variable APti + APt2 with a cut |A?7| less 
than 0.2 between the vector sums of the two pairs, 
and with 40° < |A^ < 48°. The events are gen- 
erated by the PBMl] and are charged particles of 
0.8 < Pt < 4.0 GeV/c, with |?7| < 1, from Au -f 
Au collisions at ^/SNN = 200 GeV. Since we select 
pairs of pairs which are near each other in </> (40° to 
48°) they all together pick up the bubble correlation 
and thus these pairs of pairs over all show about a 
0.4% correlation. In the -4.0 < APa + APt2 < -1.0 
region the correlation increases from 0.4% to 1%, 



while in the 1.0 < APfi +APt2 < 4.0 region the cor- 
relation decreases from 0.4% to -0.2%. This means 
the pairs of pairs are anti-aligned at a higher rate 
than aligned. This anti-alignment is what is pre- 
dicted by the color magnetic field effect presented 
above. In fact the anti-alignment increases with pt 
as the ratio bubble particle to background increases. 
However these predicted color magnetic effects have 
not been searched for yet and therefore there is no 
experimental evidence for them. Experimental in- 
vestigations of our presented four charged particle 
correlations provide a method of obtaining experi- 
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FIG. 14: Correlation lunction of pairs of pairs formed for exhibiting the effect of changes in pt due to the color 
magnetic field as defined in the text. This correlation shows that there are more anti-aligned pairs. The correlation 
is larger by ~1.2% for APti + APt2 = -4.0 compared to APn + APt2 = 4.0. APti + APt2 is equal to \ APti\ + | APt2| 
where the sign is also explained in the text. 



mental confirmation of the anti-afignment caused by 
the color magnetic field. 



IV. SUMMARY AND DISCUSSION 

In this article we have made a direct connection 
between our Parton Bubble Model (PBMpj]) and 
the Glasma Flux Tube Model (GFTM)j2]). In the 
GFTM a flux tube is formed right after the initial 
collision of the Au + Au system. This flux tube ex- 
tends over many units of pseudorapidity (Ajy). In 
the PBM this flux tube is approximated by a sum of 
partons which are distributed over this same large 
At; region. A blast wave gives the tubes near the 
surface transverse flow in the same way it gave flow 
to the bubbles in the PBM. This means the trans- 
verse momentum {pt) distribution of the flux tube is 



directly translated to the pt spectrum of the PBM. 

Initially the transverse space is filled with flux 
tubes of large longitudinal extent but small trans- 
verse size ^Qj^ . The flux tubes that are near the 
surface of the flreball get the largest radial flow and 
are emitted from the surface. As in the parton bub- 
ble model these partons shower and the higher pt 
particles escape the surface and do not interact. Qs 
is around 1 GeV/c thus the size of the flux tube is 
about l/4f initially. The flux tubes near the surface 
are initially at a radius ~5 fm. The (j) angle wedge 
of the flux tube ^ 1/20 radians or ^3°. In Sec. I 
D (also see Sec. I C) we connect the GFTM to the 
PBM (see Sec. IB): by assuming the bubbles are 
the final state of a flux tube at kinetic freezeout, 
and discussing evidence for this connection which is 
further developed in Sec. II and Sec. III. 

With the connection of the PBM to the GFTM 
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two new predictions become possible for our PBM. 
The first is related to the fact that the blast wave 
radial flow given to the flux tube depends on where 
the tube is initially in the transverse plane of the 
colliding Au + Au system. The tube gets the same 
radial boost all along its longitudinal length. This 
means that there is correlated pt among the partons 
of the bubble. In this paper we consider predictions 
we can make in regard to interesting topics and com- 
parisons with relevant data which exist by utilizing 
the parton bubble model (PBM[i|), and its features 
related to the glasma flux tube model (GFTM)[1]). 

Topic 1: The ridge is treated in Sec. II. 

Topic 2: Strong CP violation (Chern-Simons 
topological charge) is treated in Sec. III. 

We show in Sec. II that if we trigger on particles 
with 3 to 4 GeV/c and correlate this trigger particle 
with an other charged particle of greater than 1.1 
GeV/c, the PBM can produce a phenomenon very 
similar to the ridge 0, [M H [H, S M, M, m- 
(See Figs. 2-8). We then selected charged particles 
inside the ridge and predicted the correlation that 
one should observe when compared to the average 
charged particles of the central Au + Au collisions 
at = 200 GeV[ii,[3. 

In Sec. II D (Comparison to data): Triggered 
experimental angular correlations showing the ridge 
were presented at Quark Matter 2006[31|. Figure 

7 shows the experimental A(f) vs. A77 CI correla- 
tion for 0-10% central Au -I- Au collisions at y^s^N 
= 200; requiring one trigger particle pt between 3 
to 4 GeV/c and an associated particle pt above 2.0 
GeV/c. The yield is corrected for the flnite A77 pair 
acceptance. 

For the PBM generator, we then form a two 
charged particle correlation between one charged 
particle with a pt between 3.0 to 4.0 GeV/c and an- 
other charged particle whose pt is geater than 2.0 
GeV/c. These are the same trigger conditions as 
in Ref.[3l| which is shown in Fig. 7, that shows 
the corrected pair yield in the central data. Fig. 

8 shows the correlation function generated by the 
PBM which does not depend on the number of 
events analyzed. The two figures were shown to be in 
reasonable agreement when compared as explained 
previously in Sec. II D. In Fig. 9 we show the ridge 
signal predicted by the PBM for very similar data 
but with 0-5% centrality. Figure 10 shows the ex- 
traction of the jet signal. Explanations are given in 
the text. 

The second prediction is a development of a pre- 
dictive pionic measure of the strong CP Violation. 
The GFTM flux tubes are made up of longitudi- 
nal color electric and magnetic flelds which generate 
topological Chern-Simons charge [1| through the FF 
term that becomes a source of CP violation. The 



color electric field which points along the flux tube 
axis causes an up quark to be accelerated in one di- 
rection along the beam axis, while the anti-up quark 
is accelerated in the other direction. So when a pair 
of quarks and anti-quarks are formed they separate 
along the beam axis leading to a separated tt"*" and 
TT^ pair along this axis. The color magnetic field 
which also points along the flux tube axis (which 
is parallel to the beam axis) causes an up quark to 
rotate around the flux tube axis in one direction, 
while the anti-up quark rotates in the other direc- 
tion. So when a pair of quarks and anti-quarks are 
formed they will pickup or lose transverse momen- 
tum. These changes in pt will be transmitted to the 
7r+ and tt^ pairs. 

The above pionic measures of strong CP viola- 
tion are used to form correlation functions based on 
four particles composed of two pairs which are op- 
posite sign charge-particle-pairs that are paired and 
binned. These four particle correlations accumulate 
from bubble to bubble by particles that are pushed 
or pulled (by the color electric fleld) and rotated (by 
the color magnetic field) in a right or left handed di- 
rection. The longitudinal color electric field predicts 
aligned pairs in a pseudorapidity or A77 measure. 
The longitudinal color magnetic field predicts anti- 
aligned pairs in a transverse momentum or APt mea- 
sure. The observations of these correlations would 
be a strong confirmation of this theory. The much 
larger unlike-sign pairs than like-sign pairs in the 
PBM and the data; and the strong dip of the CI 
correlation at small At] (see Fig. 1 and Sec. Ill 
C for full details) shows very strong evidence sup- 
porting the color electric alignment prediction in the 
^/snn = 200 GeV central Au -I- Au collision data 
analyses at RHICfll. [l^. This highly significant dip 
(~14cr) means that like-sign pairs are removed as 
one approaches the region Ary = Acj) = 0.0. Thus 
this is very strong evidence for the predicted effect 
of the color electric field. The color magnetic anti- 
aligned pairs in the transverse momentum predic- 
tion, treated in Sec. Ill C as of now has not been 
observed or looked for. However our predicted spe- 
cific four charged particle correlations can be used to 
search for experimental evidence for the color mag- 
netic fields. 

Our success in demonstrating strong experimental 
evidence for the expected color electric field effects 
from previously published data suggests that the 
unique detailed correlations we have presented for 
searching for evidence for the predicted color mag- 
netic field effects should be urgently investigated. 
If we are lucky and the predicted color magnetic ef- 
fects can be confirmed experimentally we would have 
strong evidence for the following: 

1) CP is violated in the strong interaction in 
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isolated local space time regions where topological 
charge[l| is generated. 

2) The glasma flux tube model (GFTM) which 
was evolved from the color glass condensate (CGC) 
would be found to be consistent with a very signifi- 
cant experimental check. 

3) The parton bubble model event generator 
(PBM) is clearly closely connected to the GFTM. 
The bubble substructure strongly supported by the 
PBM is likely due to the final state of the flux tube 
at kinetic freezeout. 



ACKNOWLEDGMENTS 



This research was supported by the U.S. Depart- 
ment of Energy under Contract No. DE-AC02- 
98CH10886 and the City College of New York 
Physics Department. The authors thank William 
Love for valuable discussion and assistance in pro- 
duction of figures. 



S.J. Lindenbaum, R.S. Longacre, Eur. Phys. J. C. [19 
49, 767-782 (2007). 

A. Dumitru, F. Gelis, L. McLerran and R. Venu- 
gopalan, Nucl. Phys. A 810, 91 (2008). [20; 
D. Kharzeev, A. Krasnitz and R. Venugopalan, 
Phys. Lett. B 545, 298 (2002). [21 
L. Van Hove, Z. Phys. C 27, 135 (1985). 

S.J. Lindenbaum, R.S. Longacre, J. Phys. G 26, 937 [22' 
(2000). 

D. Kharzeev, R.D. Pisarski, Phys. Rev. D 61, [23' 
111901 (2000). 

C. Adler et al, Phys. Rev. Lett. 87, 082301 (2001). [24; 
S.J. Lindenbaum, R.S. Longacre M. Kramer, Eur. [25 
Phys. J. C. 30, 241 (2003). [26^ 
J. Adams et al, Phys. Rev. C. 71, 044906 (2005), 
S.S. Adler et al., Phys Rev. Lett. 93, 152302 (2004). [27 
X.N. Wang and M. Gyulassy Phys. Rev. D 44, 3501 [28^ 
(1991). 

K. Adcox et al, Phys. Rev. Lett. 88, 022301 (2001). [29' 
Nucl. lustrum. Meth. A 499,(2003), C. Adler et [30^ 
al., 433-436, K.H. Ackermann et al.(STAR CoL [31^ 
laboration), 624-632, M. Anderson et al., 659-678, 
F.S. Bieser et al., 766-777. [32' 
T. Sjostrand, M. van Zijil, Phys. Rev. D 36, 2019 
(1987). [33' 
J. Adams et al, Phys. Rev. C 75, 034901 (2007). 
S.J. Lindenbaum and R.S. Longacre, Phys. Rev. G [34 
78, 054904 (2008). ' [35^ 

B. L Abelev et al, arXiv:0806.0513 nucl-ex]. [36' 
J. Adams et al, Phys. Rev. Lett. 91, 172302 (2003). 

A. Adare et aZ., larXiv:0801.4020[ nucLex]. 



L. McLerran and R. Venugopalan, Phys. Rev. D 49, 
2233 (1994); Phys. Rev. D 49, 3352 (1994); Phys. 
Rev. D 50, 2225 (1994). 

T. Lappi and L. McLerran, Nucl. Phys. A 772, 200 

(2006) . 

F. Gelis and R. Venugopalan, Acta Phys. Polon. B 
37, 3253 (2006). 

S. Gavin, L. McLerran and G. Moschelli, 
larXiv:0806.4"7T8T nucLth] . 

N. Armesto, C. Salgado, U.A. Wiedemann, Phys 

Rev. Lett. 93, 242301 (2004). 

P. Romatschke, Phys. Rev. C. 75, 014901 (2007). 

E. Shuryak, Phys. Rev. C. 76, 047901 (2007). 

A. Dumitru, Y. Nara, B. Schenke, M. Strickland, 

arXiv:0710.1223hep-ph]. 

V.S. Pantuev, arXiv:0710.1882ihep-ph]. 

R. Mizukawa, T. Hirano, M. Isse, Y. Nara, 

A. Ohnishi, arXiv:0805.2795 nucl-th]. 

C.Y. Wong, arXiv:0806.2154'hep-ph]. 

R.C. Hwa, arXiv:0708. 1508 nucl-th]. 

J. Putschke, J. Phys. G: Nucl. Part. Phys. 34, S679 

(2007) . 

C. Vafa, E. Witten, Phys. Rev. Lett. 53, 535 (1984), 
Nucl. Phys. B 234, 173 (1984). 

D. Kharzeev, R.D. Pisarski, M.H.G. Tytgat, Phys. 
Rev. Lett. 81, 512 (1998). 

C. Vafa, E. Witten, Nucl. Phys. B 233, 422 (1983). 
J. Wess, and B. Zumino, Phys. Lett. 37B, 95 (1971) 
L.E. Finch et al, Phys. Rev. C 65, 014908 (2001). 



